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ABSTRACT 


This  thesis  is  a  theoretical  and  experimental  study  of  a 
scheme  for  converting  energy  from  a  streaming  fluid  through  an 
elastic  plate  to  an  electric  circuit.   The  model  is  a  quasi-one- 
dimensional  motion  of  the  elastic  plate  under  the  influence  of 
the  supersonic  fluid  and  a  transverse  magnetic  field. 

An  analysis  of  the  model  leads  to  the  equations  of  motion  of 
the  plate  in  the  unrestrained  case.   By  assuming  traveling  wave 
solutions,  the  dispersion  relation,  resulting  from  the  coupled 
equations  of  motion,  is  solved.   The  boundary  conditions  are  dis- 
cussed and  the  resulting  eigenvalue  problem  is  stated.   By  the 
method  of  coupling  modes ,  the  effect  of  the  coupling  on  the  wave 
numbers  is  determined. 

In  the  experiment,  an  elastic  media  is  tested  for  response 
with  electromagnetic  and  mechanical  excitation.   The  supersonic 
flow  is  attained  with  a  set  of  convergent- divergent  nozzles 
feeding  a  test  section.   The  combined  plate  and  flow  experiment 
is  run  and  results  discussed. 

The  conclusions  are:  (1)  the  theory  predicts  an  instability 
in  the  natural  modes;  (2)  the  plate  can  be  excited  in  the  second 
modes  and  the  frequency  response  increases  with  tension.  (3)  The 
effects  of  the  instability  in  the  combined  experiment  are  masked 
by  forced  excitation  probably  due  to  the  turbulent  boundary  layer 
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CHAPTER  I 


Introduction 

Energy  conversion  has  historically  been  accomplished  by  rotat- 
ing machinery,  making  use  of  small  gap  lengths  and  conduction  in  good 
conductors.   Recent  investigations  have  lead  to  novel  static  and 
dynamic  converters.   The  feasibility  of  generating  d.c.  power  from 
the  interaction  of  a  moving  conducting  fluid  and  an  electromagnetic 
system  has  already  been  demonstrated.     Several  schemes  for  A.C. 
power  generation  have  been  presented,  however,  these  are  restricted 
by  the  requirement  of  large  magnetic  fields  and  a  good  conducting 
fluid.   An  alternative  method  for  converting  energy  is  shown  in 
Figure  1.1.   It  is  a  device  for  producing  A.C.  power  by  abstracting 
energy  from  the  fluid  flow. 

The  device  consists  of  an  elastic  plate  in  a  channel  with  a 
fluid  flow  on  both  sides.  The  plate  is  constructed  of  highly  per- 
meable material  and  is  free  to  move  in  the  transverse  direction. 
The  channel,  also  of  highly  permeable  material,  contains  the  fluid 
and  copper  conductors.   One  set  of  conductors  is  excited  by  a  con- 
stant current  (D.C.)  source,  causing  a  transverse  magnetic  field. 
The  interaction  between  the  fluid  and  the  plate  causes  transverse 
motion  of  the  plate,  similar  to  the  "flutter"  motion  in  aero-elasti- 
city. This  motion  provides  a  variable  coupling  of  the  first  cir- 
cuit with  the  second  magnetic  circuit,  producing  A.C.  power,  as 
indicated  in  the  lumped  parameter  model  shown  in  Figure  1.2. 

Analytically  the  problem  of  the  thin  airfoil  performing  small 
lateral  oscillations  in  a  uniform  stream  of  incompressible  fluid 
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was  for  many  years  the  heart  of  all  flutter  prediction.   The  com- 
plete solution  was  first  published  in  the  United  States  by  Theo- 

(2) 
dorsen   .   This  approach  has  been  further  developed  by  many  workers 

and  is  presented  most  comprehensively  in  the  text,  Aeroelasticity . 

(3) 
by  Bisolinghoff ,  etc. 

The  definition  of  panel  flutter  is  the  self  excited  oscilla- 
tion of  the  external  skin  of  a  flight  vehicle  when  exposed  to  an 
airflow  on  one  side   The  analysis  of  panel  flutter  has  been 
pursued  in  great  detail  recently,  but  there  is  still  some  disagree- 
ment between  theory  and  experiment.   It  appears  that  the  theoreti- 
cal considerations  to  date  do  not  approach  the  situation  described 
here.   The  tensions  on  the  panel  are  significantly  less  than  those 
to  be  applied  in  this  device.   The  mechanism  for  flutter  is  not 
explained  physically  in  the  context  of  a  feedback  situation.   How- 
ever, the  effects  of  aspect  ratio,  elastic  foundations,  damping 
and  boundary  support  have  been  covered  in  detail.   In  the  appendix 
a  comparision  of  the  prediction  of  flutter  theory  and  the  pre- 
dictions of  this  theory  is  shown.   The  aeroelastician  in  general 
has  confined  his  analysis  to  that  of  "infinite"  flow,  whereas  in 
this  scheme  the  confining  effect  of  the  rigid  walls  will  be  impor- 
tant.  In  dense  fluids,  there  has  also  been  investigation  of  flow 
induced  vibrations  of  fla 
followc 
recent . 

elastic.   *    The  latter  is  concerned  with  a  statistical  repre- 
sentation of  the  turbulent  boundary  layer  as  a  forcing  function. 
Here  a  rain,  the  models  have  infinite  geometry  and,  therefore,  are 
not  directly  applicable   Lastly  there  are  several  "handbook"  solu- 
tions available  for  determining  natural  modes  of  vibration,  for 
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example ,  see  Den  Hertog    . 

(8) 
Melcher    has  studied  field  coupled  surface  waves  for  both 

electrohydrodynamtc  and  magnetohydrodynaroic  systems.   Ketterer, 

studying  under  him,  investigated  the  interactions  between  a  fluid 

stream  and  a  fixed  structure  capable  of  supporting  traveling 

waves   By  theory  and  experiment,  it  was  shown  that  energy  could 

be  extracted  from  the  flow  in  the  form  of  oscillations. 

In  this  study  there  is  contact  between  the  fluid  and  the 
elastic  structure  which  was  not  present  in  the  experiment  conducted 
by  Ketterer.   The  effect  of  this  contact,  a  boundary  layer  which 
could  be  turbulent,  may  be  considerable  in  comparison  to   the 
effect  of  an  instability  in  the  natural  modes. 

There  are  several  methods  of  exciting  the  plate  available. 
First,  it  could  be  located  transversely  in  the  flow  path,  simi- 
lar to  holding  a  piece  of  grass  between  the  thumbs  and  blowing. 
However,  the  length  of  interaction  would  be  limited  and  therefore, 
the  amount  of  energy  taken  from  the  flow  would  be  limited. 
Secondly,  the  plate  could  be  fixed  at  one  end  only,  and  the  exci- 
tation would  be  dependent  on  the  trailing  vorticies.   This  would 

(4) 
be  similar  to  the  study  reported  by  Eagleson.     Another  method 

that  has  been  suggested  is  to  place  the  plate  in  a  tube  closed  at 
the  downstream  end  and  put  the  combination  in  the  supersonic  flow. 
This  would  provide  greater  excitation  since  pressure  variations 
considerably  above  stagnation  pressure  have  been  observed  in  such 
a  configuration.   The  overall  efficiency  of  this  scheme  is  ques- 
tionable, however. 

The  method  chosen  was  considered  both  theoretically  and  experi- 


9H 

taa»i 
■0V0M 

■ 

it1     :  but*    «i< 

)fl  &•  for  gftu  i 

•    MIX 


:>*na&!a    9vr   ' 


. 


I 


mentally  the  simplest  to  study,  with  the  added  feature  that  the 
significant  interaction  length  might  produce  a  reasonably  effi- 
cient devLce.   In  Chaoter  IT,  the  device  is  studied  by  assuming 
a  quasi-one  dimensional  motion  of  the  plate.   The  resulting  dis- 
persion relationship  is  solved     aaoter  I  EI  and  the  inolica- 
tions  discussed.   In  chaoter  IV  the  longitudinal  boundary  condi- 
tions are  treated.   The  experiment  and  results  are  described  in 
Chaoter 
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CHAPTER  II 


ANALYSIS  OF  DEVICE 


The  device  will  be  modeled  by  assuming  a  quasi-one  dimen- 

sional  motion  of  the  plate.   The  model,  as  seen  in  cross -section 

from  a  top  view,  is  shown  in  Figure  2.1.   The  flow  will  be  assumed 

inviscid  and  incompressible  in  this  analyses  since  the  objective^ 

I 
is  a  device  that  can  be  tested  in  a  wind-tunnel.   The  proposed 

orientation  will  also  allow  neglecting  the  effects  of  gravity. 

The  fluid  flow,  the  nlate  and  the  magnetic  effects  will  be  treated 

separately,  and.  then  the  results  wil     combined  to  obtain  the 

equations  of  motion  ^ST 

2  1  lie       Jquatlons! 

In  the  bulk  of  the  flg^U  ^  body  forces  are  neglected,  the 
momentum  equation  is 
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4na  lytical 


for  an  inviscid  fluid  and  the  equation  for  continuity  of  mass  is 

^  +  p  V-y.  =  0  (2.2) 

Dt   r    l 

FKkURB    1.1 

for  compressible  flow    '.*  will  assume  small  perturbations  and 

that  the  waves  are  long  compared  to  the  significant  dimensions, 
that  is: 
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v**   -  v'  (x,t)  -  V 
x    x   *      o 


v'   ■  v  (x,t) 

y   y 


p(x,t)  -  p'  (x,t)  -  pq  <  <  oq 


p(x,t)  -  p'(x,t)  -  P0  <  <  PQ 


where  the  primed  quantities  are  the  actual  values  and  the  unpri- 
med  quantities  the  oer tar bat ions     serting  these  into  equations 
(2  1)  and  (2.2)  arid  retailing  oql^  finest  order  terms: 


^v 
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Initially,  the  thermodynamic  process  will  be  considered  adia- 
batic,  which  allows  relating  the  ebanges  in  oressure  and  density 

c  ciRcur 

by  a  constant. ____——. 
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2.2  Magnetic  Equations 

A  section  of  the  magnetic  circuit  is  shown  in  Figure  2.2. 

We  assume  a  perfectly  magnetic  material  so  that  the  coil  will 

induce  a  magnetic  field  in  the  circuit,  with  the  H  field  existing 
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in  the  gap  only.   From  Ampere's  circuital  law  we  get: 
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Define:   H  *  -r-  as  the  steady  state  H  field.   Therefore,  the 
o   a 

magnetic  fields  with  reference  to  figure  (2.1)  will  be: 

H  '  -  H   (  I — - —  ) 
y 


where  the  superscripts  refer  to  the  resoective  sides  of  the  slate 
There  are  several  methods  available  for  computing  the  force  on 
the  plate.   The  most  direct  method  is  to  use  the  Maxwell  Stress 
Tensor  described  in  detail  in  reference.     The  force  density 
in  the  y  direction  acting  on  the  plate  is  the  sum  of  the  stresses 
acting  in  that  direction  on  each  of  the  faces  of  the  plate.   The 


stress  tensor: 

as 


In  our  case,  where  H   is  the  only  significant  component, this 
becomes 
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and  the  force  per  unit  area  is 
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Keening  only  the  linear  terms,  the  irorce  density  is  simoly: 


Therefore,  the  force  is  proportional  to  the  displacement  and  will 
have  a  destabilizing  effect  on  the  plate. 

2  3  The  Plate  Squat ion 

The  plate  has  a  density  p  and  a  modulus  of  elasticity  E. 
Using  simple  bea<6  theory,  the  bending  moment  at  any  point  is  pro- 
portional to  the  rate  of  change  of  the  slope. 

2  3 

m  ■ 1     and  dm  = if   dx 

w  ^  2  w  ^  3 

^x  dx 


Taking  an  elemental  section  as   in  Figure2.3    we  can  indicate  all  of 

the   forces  and  moments  acting  on  it. 
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f  ■  external   force   density   =  2jjl  H     ■*  4-   (p   -p   ) 

J  o  o  d    r  r 

v  *  shear  force  per  init  width. 

p  -  tensile  force  per  unit  width. 
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m  ■  bending  moment. 

Now,  summing  the  moments  about  the  lower  left  hand  corner: 
-dm  -  vdx  +  pdi  «=  0     or 


dm  ,    Ji      EI  d3    ,    il 
dx     dx      w   .  3    v   dx 


and  the  forces  acting  on  the  element. 


-   T~  ~  pmzA  — *   or 
at 

2  4^2 

f  +  P  — o  -  —  d  — f   -  pm2A  — * 


w    dx4        at2  (2.5) 


The  three  previous  developments  can  be  combined  to  give  the  equa 
tions  of  motion. 

2.4  The  Equations  of  Motion 

Equation  (2.4)  can  be  rearranged  and  integrated  directly. 
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The  transverse  fluid  velocity  must  be  zero  at  the  channel  wall 
in  both  the  uooer  and  lower  regions: 
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(2.7) 


If  one  assumes  no  cavitation  the  transverse  fluid  velocity  must 
equal  the  plate  velocity  at  the  *>late. 
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the  above  two  equations  maybe  added: 
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Equation  (2.7),  the  plate  equation,  maybe  rearranged  and 

u  ■ 
inserted  in  (2  3)  and  (2.8)  to  eliminate  (p  -p  ),  the  pressure 

difference,  from  our  equations*. 
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x  ox  >t 


(2.9) 


,  i.  +  v     |-  )t  -  — *-  (  ^-  +  v     i-)(2aH2-.^^  +  o^ 
<   >t  +  vo  "x   )?        0        2    (   ot       Vo  SxMV     d       w     .   4  +  P  ox 

2p   c  ox 

-   2pmA  -— |  )    "   2  ^x     v  =  ° 

(2.1C) 


where  we  have  defined  V  ■   (V     -  V  ) . 

x         x 
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It  will  be  convenient  to  have  the  equations  in  nondimensional 
form  before  further  manipulations : 

Variables:   n  i¥  Z/d,   v*  ^-  ,  x  *  x/d,  k  -  kd,  <L  •  a>  -i-  ,  t  «=  ■—  t 

v  v  d 

o  o 

P  d 
o 
Parameters:  n   *  — r  , density  ratio, 

M     =   ~    ,Fluid  "Mach"  number 

a        v      • 
o 

1/2 
M     .   (P/2Apm)    '       ,  Material  "Hach"   number 
m  v 

o 

HoH~d 

H  ■  r-   ,  Magnetic  pressure 

pmvQA 

d  =  u — 2T"  3  Risidity 

w2prov  d  A 
0 

We  have  normalised  the  parameters  to  a  dynamic  pressure 

1    2 
(  -r  p  V  >  and  then  multiplied  by  the  density  ratio.   This  will 

gimol...  fy  the  solutions  as  well  as  provide  an  indication  of  the 

strength  of  the  various  forces  or  pressures  compared  to  the  fluid 

dynamic  pressure.  « 

dpoc 
We  now  multiply  equation  (2.9)  by  r  and  equation  (2.10) 

ap  \r 

rm  o 


T: 


i 


VfD< 


>  3  »mS7  »S 


c 


by 


d2 

2Ap  v2 
rm  o 


nxf  <  Is  +  |s  )  n-(  Is .+  Is)  <H-„  M    +  m2  2l_  v    *  „2  1=. 


(2.11) 


(2.12) 


The  above  form  the  equations  of  motion  of  the  system  in  nondimen- 
sional  form 

To  solve  these  equations  as  they  stand  for  specified  forcing 
function  would  be  difficult  indeed.   This  approach  has  been  fol- 
lowed on  simplified  equations,  where  the  forcing  function  was  a 
statistical  representation  of  the  boundary  layer  excitation.   The 
approach  we  will  follow  is  to  look  for  instabilities  in  the  free 
modes,  when  the  source  of  oscillations  is  not  a  lossy  one  as  it 
is  in  a  turbulent  boundary  layer. 
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CHAPTER  III 


THE  DISPERSION  EQUATION 

i 


The  dispersion  equation,  or  the  relationship  between  d>  and  k, 
will  give  an  indication  of  the  stability  of  the  plate  in  the  free 
modes.   We  will  assume  a  delta  function  in  time  and  space  as  the 
excitation,  so  that  the  equation  of  motion  will  have  no  driving 
function   In  order  to  determine  whether  the  structure  can  sup- 
port traveling  waves,  we  also  assume  that  the  solutions  will  be 
of  the  form^  ^  * '  e3  (a)t-kx)        andv  m  ^   eJ(a>t-kx).   Insert- 
ing these  into  the  equations  of  motion  we  get: 


(i-k)(ntf2  -  H+DkViV^2)^  +  ^  (K2k)v  -  0 

3  (3.1) 


-k  (H-Dk4-  I  2k  +  J)2)^  +  ^  (i-k)^  -  0 


m 

(3.2) 

Setting  the  determinant  of  the  coefficients  equal  to  zero 
will  insure  that  the  equations  of  motion  are  compatible: 


ir  {  (o>-k)?(uK2-H+Dk4  +  K2k2-  m2)  +  M2k2 (H-Dfc4-M2k2-f  a>)}  =  0,„  „. 
2  w  *■  a  m  a         m      ' 
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This  is  the  desired  dispersion  relationship,  and  its  solu- 
tion, 05  *  f  (k)  ,  and  the  inverse,  k  -  f  («e)  ,  will  provide  insight 
into  the  stability  of  the  infinite  plate.   Before  proceeding  with 
the  solution,  there  are  several  limiting  cases  which  should  be 
investigated.   These  in  effect  will  separate  the  two  systems. 

1)   If  we  allow  the  tension,  P,  on  the  plate  to  become 
infinite,  the  plate  becomes  a  rigid  boundary. 
Thus,  if  M  — *o»  ,  we  get: 

.  -  9-9     9-9    -9 
(cjo-k)V  +  MTc   (-k^)  -  0 

or,  in  dimensional  form: 


(<d-vow2  =  mV  (Ji4) 


which  is  the  expected  dispersion  equation  for  convective  acoustic 
waves 

2)   Similarly  if  we  allow  the  density  of  the  fluid  to  approach 
zero,  we  find  the  equation  of  the  plate  alone. 


.■>  -4    2-2  <35> 

a)-  -  -  H  +  Dk  +  M  k 

m 
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&  aril   ai   t 
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It  is  comforting  to  note  that  the  magnetic  force  will  have  a 
destabilizing  effect,  while  the  tension  and  rigidity  tend  to  sta 
bilize  the  solution. 

3)   Allowing  the  speed  of  sound  to  become  infinite  will  res' 
trict  the  solution  to  the  subsonic  region.   In  dimensional  form: 


2u.H 


I   w  ^  4       2  ~  2    ho  o  d 

^K        -'X  <^t 


This  is  the  equation  for  the  plate  with  a  dynamic  pressure  term 
which  can  be  derived  from  Bernoulli's  equation  in  a  linearized 
form  (as s>iming  tt  — >> 

P   4  ~  r  V?  =>  ?  4 \a{   -r™r  )   -"'?   i-  P'  +7  pV   4-  .pV2  h 
o   2  *  o  o       2^ v  d-^       o        2  H  o»    h  o»d 


then: 

Pu-P  '  -  2P  V2  | 
ro  o   d 

The  above  three  checks  indicate  that  the  model  and  the  ana- 
lysis are  consistant.   It  also  appears  that  they  are  to  some 
extent  physical. 

The  third  check  gives  a  static  instability  which  has  been 

(10) 
observed  at  subsonic  speeds. 

The  uncoupled  equations  (3.4)  and  (3.5)  can  be  solved  simply, 
since  they  are  quadratic,  and  plotted,  as  in  Figure  3.1.   For  the 
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rest  of  the  analysis,  some  experimental  considerations  will 
govern  the  approach.   Supersonic  soeeds  are  difficult  to  attain, 
and  in  general,  the  experimental  apparatus  is  designed  for  a 
specified  flow   Therefore,  we  will  take  V  to  be  a  constant 
equal  to  twice  the  speed  of  sound.   The  rigidity  of  the  plate  is 
also  a  fixed  quantity  for  a  given  experiment  and  the  effect  of 
its  variation  would  be  roughly  equivalent  to  varying  the  tension 
on  the  plate.   Since  the  plate  tension  is  easily  changed  experi- 
mentally, it  will  be  allowed  to  vary  in  the  solution  of  the  dis- 
persion equation. 

We  now  rearrange  the  dispersion  equation  as  a  oolynomial  in 

0)  •. 


^4-to3(-2kH  u2   (H-Dk4-  M2k2-  M2-k*M2  +  k2  ) 

m     a    a 

+  &   (!>k)  (-H  +  Dk4  +  M2k2+  M2) 

III  O 


-9    9  -L  9-9 

+  (-k^)(HZ  -  1)(H-Dk*  -  MHc  )    -  0 
a  m 


(3.6) 


and,  using  the  computer,  solve  for  real  values  of  k 
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The  results  are  shown  In  Figures  3.2,  3.3,  and  3.4  for  the 
following  cases: 


1)  Figure  3.2 

2)  Figure  3.3 

3)  Figure  3.4 


M   -  .875 

m 


m 


545 


M   -  .314 
m 


(strong  tension) 
(medium  tension) 
(weak  tension) 


The  solutions  plotted  in  Figure  3.1  remain  as  the  asympotes 
for  these  solutions  Returning  to  equations  3.3  and  3.4,  we  find 
in  dimensional  form: 


cjd  f  (7±c)k 
1/2 


for  th 

le   fluii 

05=^  ( 


for  the  olate 


In  the  case  of  strong  tension  the  coupling  is  weak.   Kesrever,  as 
the  tension  is  decreased  as  in  Figure  3.3  to  a  point  where  (v  ,), 
the  wave  speed  in  the  plate,  is  approximately  equal  to  (V-c) , 
the  "slow*  fluid  wave,  the  coupling  becomes  pronounced.   Decreas- 
ing the  tension  further,  decreases  the  coupling  until  the  solu- 
tions are  uncoupled 

The  two  systems,  the  fluid  and  the  plate,  are  in  themselves 
stable  at  at  least  if  the  majne  tic  field  is  small.   The  above 
solutions  indicate  that  when  coupled  together  an  unstable  situa- 
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tion  is  produced.   That  is:   cd  has  a  negative  imaginary  part 

ir   „  ,  \  it+j  (a>t-kx) 

and  the  assumed  solution:  n==ll  eJ  V.a>t-Kx;   has  the  form:  n  c  ne 

such  that  the  displacement  grows  in  time.   This  is  the  result 
of  a  feedback  mechanism  within  the  coupled  system.   Let's  assume 
a  deflection  of  the  plate  in  the  second  mode,  with  a  positive 
amplitude  over  the  leading  half  of  the  plate.   This  provides  a 
compression  wave  in  the  fluid.  As  this  disturbance  is  swept 
downstream,  it  arrives  over  the  trailing  half  of  the  olate  in 
time  to  reinforce  the  motion  of  plate.   The*-*  the  motion  is  trans- 
mitted 


-t  --  o  L  s  4 

back  upstream  through  the  olate,  to  arrive  at  the  right  time  to 
start  the  orocedure  over  again.   This  is  the  case  when  the  velo- 
cities /  .  and  V-c  are  equal.   Since  this  is  a  distributed  oheno- 
menon  over  the  half  length  of  the  olate,  the  mechanism  will  work 
for  velocities  that  are  approximately  equal  as  well  as  exactly 
equal. 

A  complete  discussion  of  the  dispersion  equation  as  it  per- 
tains to  the  infinite  problem  would  not  add  significantly  to  the 
study  of  the  device  proposed.   However,  the  steps  will  be  outlined 
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Having  solved  for  a)  *  f(k   ..),  the  complex  values  of  oi  obtained 

real 

are  substituted  and  the  equation  solved  for  k  +  jk,  ■  f(o>  +  jo),) 
to  determine  the  type  of  instability  involved.   This  can  be  done 
using  the  Bers-Briggs  criteria,     with  the  help  of  computer 


technique  givi.ag  a  graphical  display  of  both  co  versus  oj.  and 

(12)       r 
k  versus  k.  as  described  by  hills. 
r        1  7 

The  disoersion  eouatlons  shown  in  Figures  3,2  -  3.4  have 
characteristics  which  indicate  that  the  basic  theoretical  model 
is  physically  meaningful   The  coupling  between  the  flow  and  the 
elastic  media  occurs  at  long  wave  lengths,  with  the  short  wave- 
lengths stable  and  uncoupled,   These  are  the  ^ronerties  of  the 

stream  structure  previously  investigated  in  the  context  of  an 

(8) 
electrohydrodynamic  interaction      There  it  was  shown  that  the 

effect  of  the  longitudinal  boundaries  was  to  make  the  instability 
absolute.  That  is,  the  theory  supports  the  conjecture  that  grow- 
ing oscillations  can  be  oroduced. 


ft 


■     91 B 


.1,1 


ai 

1;     8*7  Slbftfll    ollB&I* 

';   9B*  n»  sfdelE    erto$n9J 

ICf     9"  i9T38 

JO»J» 

09^5» 

■i    1r  9"JrjJ  08<J» 


31 
30 


CHAPTER  IV 


BOUNDARY  CONDITIONS 


The  a?*)  I  Lett  ion  of  the  boundary  conditions  will  be  aporoachei 
by  two  methods,  an  exact  method  resulting  in  the  eigenvalue  pro- 
blem and  an  approximate  method  called  the  coupling  of  modes 
approach.   A  sitinle  linearized  example  will  be  covered  first  fol- 
lowed by  a  discussion  of  how  to  apply  the  boundary  conditions. 
This  leads  to  the  statement  of  the  eigenvalue  problem.   Next,  the 
coupling  of  modes  is^4ttf«[I&»*4^  This  will  give  an  indication  of 
a  more  meaningful  wmj  ut   mrtytng  the  boundary  conlitl'dffa',  and 
an  estimate  of  the  coupling  effect.   The  boundary  value  problem 

is  not  solved  completely. 

TENSILE     STRESS  -  P 
DENSITY  -  0 

4  1    A  Simple  Example 

Consider  a  membrane  fixed  at  the  ends.     This  is  similar  to 

the  elastic  plate  if  the  rigidity  and  the  fluid  are  neglected. 
A   MEMBRANE    F/tfEO     AT   Trtfe   ENDS 


2  2 

I       Equations   of  motion:  d   I*  d  £ 

P  ftt2   "  P  dx2 

2,     Assume  Solutions   of  the   form:        $   ■   (A  sin  k     +3  cos  k  )   e^ 
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A   MEMBRANE    FIXED     AT   Ttt£   ENDS 


FIGURE     4.1 
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4        Boundary  Conditions: 

l)  ft*,  A  - 


a 


2)      ^(Ltt)    *    j       A  sin  kL  «     and   k  L  331 
Let 
5.   Solution:  1 


5(x,t)  -Re   f  (An  sin  *f  )  e^  ((  f,     f  )  «=  J 


4  2 -  Method  of  Characteristics 

In  the  example  the  basic  steps  to  the  oroblem  solution  are 
stated.   There  is  a  significant  change  when  the  fluid  is  aided, 
especially  when  it  is  allowed  to  have  a  grog's  translational 
motion   The  question  arises  -at   what  noint  may  one  at>p  ly  boun- 
dary conditions  on  the  flow  such  that  they  are  meaningful?   Jere 
the  fluid  motionless,  the  boundary  conditions  must  be  applied 
at  two  ooints  in  st>aee.   This  is  also  true  if  the  velocity  of  the 
fluid  is  less  than  the  speed  of  sound  in  the  fluid.   From  the  dis- 
persion equation,  however,  the  instabilities  occur  at  sxipersonic 
speeds  only. 

In  order  to  understand  this  more  clearly,  the  method  of 
characteristics  is  useful.   This  method  uses  the  general  solutions 
of  the  wave  equation  for  the  transient  case: 

f  -  f  +  (x-cO  +  f  -  (x+  ct) 
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where  ff  is  a  waya  propagating  in  fha  p1ni,x  direction  and  f-  a 
wave  in  the  negative  x  direction.   The  solutions  are  represented 
by  thf  third  dimension  of  a  $raoh  where  the  first  two  are  time  and 
space.   Thus  a  disturbance  can       '  <STT-|lJr both  time  and  space, 

Let  :    i  =  x  +  ct   and  B  -  x  -  ct 

FIGURE   4-2 

Assume  a  disturbance:      P    -  a  <   <   a 

o 

•  0  elsewhere. 

for  the  fluid  flow  shown  in  Figure  4.,^. 


L. 


In  Figure  4.3  the  disturbance  is  showriT* 

! 

Lines  representing  constant  qp  and  B  can  be  drawn  on  the 
graph   Those  originating  at  x  ■  +  a  are  the  most  interesting 
since  they  show  the  progress  of.   the  disturbance.   For  example, 
at  t  =  t  the  position  of  the -disturbance  cart  be  determined  from: 

f  -  f+o>    I-  t0  +  e-/«)     t  -  to     /  ; 

The  tlooe  of  the  lines  is  of  interest  since,  if  the  lines  of  B 
had  a  positive  slooe,  the  point  x ■  •=  -L,  would  never  know  the  dis- 
turbance  This  is  the  case  when  the  fluid  is  allowed  to  have 
motion  at  a  velocity  greater  than  the  speed  of  sound.   For  a  trans 

lating  media,  we  must  redefine: 

FLUID    WAV/E  PROPACAT 

i  ■  x  +(v+c)t  and  8  =  x+(v-c)t 

FIGURE     4.3 
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as  the  characteristic  solutions  to  the  convecting  wave  equation. 
Then  (v-c) ,  the  slope  of  the  9  lines  becomes  positive  for  v  >  c, 
and  points  upstream  will  not  be  effected  by  events  downstream. 

This  implies,  in  the  scheme  proposed,  that  the  fluid  boun- 
dary  conditions  must  be  applied  at  the  leading  edge  of  the  ^>late 
in  order  to  determine  their  effect  on  the  interaction.   The  fluid 
will  not  know  what  is  happening  at  the  trailing  edge  until  it 
gets  there.   There  are  two  variables  on  the  flow,  the  pressure 
and  velocity,  and  their  values  at  the  ooint  x  ■  0,  as  a  function 
of  time,  form  the  constraints  on  the  flow. 

4.  3  The  Eigenvalue  Problem 

The  dispersion  equations  was  sixth  order  in  to   .      The  six 
boundary  conditions  are  the  two  specified  for  the  flow,  plus  the 
four  on  the  plate   From  simple  beam  theory  the  latter  are  deter- 
mined from  the  type  of  support-free,  pinned,  or  clamped.   If  we 

assume  the  ends  are  clamped,  this  implies  that  that  the  deflec- 

d£ 
tion  and  slooe,  -p  ,  are  zero  at  x  ■  0  and  x  -  L.   Since  v(x,t) 

and  p(x,t)  are  defined  as  difference  quantities  (above  and  below 

the  plate"),  their  values  are  zero  at  x  =  0. 


Returning  to  the  equations  of  motion  this  can  be  stated  sim- 
ply. 

(©  -  k)(|iMa-?.E.)r  •:   |  (Ma?k2)  v  -  0  (4.1) 

k(P.E.)"  4-  |  (m-k)  v  -  0  (4.2) 
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(P.E.)r)  +0  p  -  0  t  (4.3) 


where  P.  E.  -  H-Dk4-M2k2  +  a)2 

m 


As  in  the  example,  assume  a  solution  of  the  form: 
r,(x,t)  -Re  (J  Ae^V^} 

and  rewrite  equation  (4.1)  and  (4.2). 

P  -  -  ~     D  (4.4) 

v  -*-  Xffi=£i—  t,  -  £f4  (4.5) 

1        k(P.E.) 

apolying  the  boundary  conditions: 

(1)  T)(0,t)    - .0  i  Ajj   =  (4.6a) 

(2)  |«   (0,t)    -  0  ^  JHi^    =  0  (4.6b) 

At 

(3)  n-(L,t)    =0  £    A^   ejknL  -  0  (4.6c) 

(4)  £-     ri(L,t)        -0      £2ik  %    ejknL  N  0  (4.6d) 


(5)      p(0,t>    =0  £  2(   -    (p>E]-(ktr)   *rf   -  ^  (4.6e) 
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(6)   v(0,t)  =0       £  2f  (k  )  A  «  0  (4.6f) 

L      n        n 


Rearranging  in  matrix  form: 
A(J,k)~ 


1 

1 

1 

i 

1 

1 

Ri 

k2 

S 

H 

S                    ^6 

%l 

«fi 

e^3L 

e^4L 

e-'    3 

e^6L 

££■£ 

k2a^2L 

k3eJ£3L 

k.eJ   4 
4 

k5eJSL 

,-      ik,L 
k6e     6 

1 

1 

1 

I 

1 

1 

P.E    (£   ) 

P.E.(k"2) 

P.E.(k"  ) 
3 

P.E.  (If.) 

P.E.(k_5) 

P.E.(k6) 

ft(kt) 

,(k2) 

fi(V 

vv 

£i<V 

f,(k6) 

In  theory  then  the  equations  ^(a}k)  =  0  and  D(co,k)  =  0  can 
be  solved.   Assuming  a  complex  cfi,  the  six  k's  are  found  from 
D(cu,k)  ■  0.   These  values  are  substituted  in  to  A.   This  procedure 
is  continued  until  A  =0;  both  it  s  real  and  imaginary  parts, 
since  it  is  complex   In  practice  another  method,  that  of  weak 
coupling  of  modes  will  be  used  to  attack  the  problem. 
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4.4  Method  of  Coupling  of  Modes 

One  might  ask  the  question:   'Tf  two  systems  are  individually 
stable,  why  would  a  weak  coupling  between  the  two  result  in  totally 
different  motions?'   The  answer  is  that  the  basic  solutions  for 
the  two  systems  are  approximately  correct,  and  that  the  changes  can 
be  predicted  within  certain  limitations.   First,  the  systems  must 
be  stable  by  themselves.   In  this  problem,  this  restriction  will 
not  allov   including  the  effect  of  the  magnetic  stress  on  the 
plate  since  it  has  destabilizing  influence.   Secondly,  the  coup- 
ling between  the  two  systems  must  be  weak.   Within  these  restric- 
tions,  mutual  coupling  of  modes  will  allow  simplifications  based 
on  a  physical  understanding  of  the  system. 

(13) 
In  theory,  the  approach  is  quite  simple.      A  complicated 

system  is  first  divided  into  its  component  parts.   The  equations 
of  motion  for  the  isolated  systems  are  solved  exactly,  and  the 
solutions  are  expressed  in  the  'normal  modes   of  the  element. 
The  complex  system  is  then  assriraed  to  be  made  up  of  these  ele- 
ments weakly  coupled.   The  motions  within  the  system  are  described 
by  the  elemental  motions  which  are  slightly  perturbed  by  the 
coupling. 

As  an  example,  in  the  system  under  consideration  the  ele- 
mental parts  are  the  fluid  and  the  plate.   Assuming  the  fluid 
motionless  and  the  plate  of  negligible  rigidity  (a  membrane) , 
each  will  propagate  disturbances  as  waves,  both  a  forward  wave 
and  a  backward  wave  in  each  system.   As  seen  from  the  solution 
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of  the  dispersion  relationship,  the  coupling  will  occur  between 
the  forward  membrane  wave  and  the  "backward'  fluid  wave.   There- 
fore, the  other  two  waves  can  be  neglected.   The  motion  of  the 
system  will  be  that  of  these  two  waves  slightly  altered  by  the 
coupling. 

4. 5  Elemental  Equations  in  the  Coupled  Mode  Form 

JL.   Fluid  Equations 

The  fluid  equations  can  be  greatly  simplified  to  get  an 
approximation  to  the  solution. 

^v(x.t)  _      1     ^o(x,t) 

A.  Zx  P^(v+c)   dt  (4,8a) 

o  ~~ 


NP<*,t)  .     Mx,t) 

^x         po  *t  (4.3b) 


These  coupled  equations  are  most  easily  solved  if  a  method  of 

(13) 
decoupling  them  can  be  found   Following,  Louisell,     the 

decoupled  form  is  derived  from  linear  combinations  of  the  above 

equations   That  is;  assume  a  normal  mode  form: 

1 


b  +  <x,t>  -  \  (P)  (V(x,t)  +  p  (^c)  P(x,t)) 
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where  the  mode  amplitudes  are  normalized  so  that  the  square  of 
the  amplitudes  represents  energy  carried  by  the  mode. 

Multiplying  equation  (4  3b)  by  the  fluid  impedance  (  — ,    .  ■>  ) 

po  ~ 
and  adding  and  subtracting  from  equation  (4.8a)  will  give  the 

desired  decoupled  form: 


x  ^  ±  P  (V+c) 


p)  -  -    ?fe'tt**  infer  p) 


or 


^  (b±  (x,t))  «  -  -tLI:  (b±  (x,r)  ) 


^x 


V+c  -)t 


Assume  that  the  time  dependence  will  be  of  the  form  eJ   , 
then:  b  +  (x,t)  -  a  +  (x)  eja)t  +  a*+(x)  e~ja)t 


Now  if  v(x,t)  -  Re  V(x)e:5a3t 

1 

then  a  +  (x)  -  |  (p)2(V(x)  +   (v+c)   P(x)) 

o  ~~ 

and  equations  (4.8a)  and  (4.8b)  become 


(  ~  +  jk)a+(x)  -  0 


(  $T  "  jk)a*  (x)  *  0 


(4.9) 


(  -r-   -  jk)a-(x)  -  0 


r:~-+  1k)a*(x) 


(4.9b) 
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where     k  =  — ; — 

vH  c 


2       The  P7ate  Equations 

Again  simplifying  assumptions  must  be  made.      Neglect- 
ing the  rigidity  and  the  external   forces  on  the  plate,    the  nor- 
mal mode   form  of  the  equations    for  this  element  can  be  derived 
in  an  analogous  manner. 


^x  p  ot       *t  Wtiere  vu'c;        }t 

*t   (Xjt)        p   ^x  bU,t>        *x 


then  i 


b  +    (x,t)    -  \  (p)7(  7(x,t)   +  §  S(x,t)) 


and  i 

2 


fltj     (b±  (x,t))    *%($f£     (b±  (x,t)> 


there  fore , 

(  r^  -jk2)a2+(x)   =0  (  |j  -jk2)   a*  +  -  0  (4.10a) 

(  |r  +jk2)   «2"(k)   -  0  J:|j  -jk2)a*  -  -  0  (4.10b) 
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where 


k,  --a- 
(p/p) 


4- 6  Coupling  Between  the  Two  Elements 

In  section  B,  we  derived  the  normal  mode  form  of  the  element 
equations    In  effect  the  single  variables  in  equations  (4.3a)  and 
(4  3b)  have  been  replaced  by  a  combination  of  the  two  variables. 
To  put  the  previous  equations  in  coupled  mode  form,  it  is  nec- 
essary to  make  linear  combinations  of  the  elemental  forms.   This 
is  done  in  such  a  manner  that  when  the  couoling  is  reduced  to 
zero,  the  equations  reduce  to  those  derived  for  each  element 
separately. 

'^al 

^T  "  cllal  +  C12a2  +  c13al  *  C14a2  (4'lla) 


^a 


IT   =  c21al  4   c22a2  *   c23al  +  C24a!  (4.11b) 


I 

5*         ^31al    '    "32a2     "33aI    '    "34a2 


c^a,  +  c~na  +  c^„a*  +   c0/a*  (4.11c) 


>■-!* 


■  c,.,a,  +  c/na~   +  c/0a?  ,  c.,as  (4. lid) 


^x    ^41° 1  '   42a2   "43~1  +  "44"2 


where 


'i-j(^V>  +  F"fe)   ?(x)> 
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l2  ■2(V/tV?)  fpS<»») 


and 

c22  =  -jk2(l  +  t/4Ffk2) 


c   -  -1  -J— 
°32    J  4pfk2 


o   A  coupling  coefficient 

In  order  to  make  the  problem  more  tractable  we  must  make  use  of 

our  previous  results.   The  first  is  that  coupling  between  the  two 

systems  will  occur  when  the  two  wave  members,  k,  and  k« ,  are 

approximately  equal.   Since  k-  ■  — ^—  and  k«-7r-  and  in  Chapter  II 

m 

we  saw  that  the  coupling  occured  when  /  +  =  V-c,  we  can  there- 

m 

fore  neglect  the  fast  fluid  wave  (V-tc)  and  the  backward  plate 

wave  (V  ) . 
m- 

Secondly,  we  must  determine  if,  in  fact,  the  coupling  between 
the  two  systems  is  weak.   Returning  to  the  coupled  dispersion 

relationship  and  inserting  our  previous  assumption  that  the 

rigidity  and  the  magnetic  field  are  zero,  we  have: 
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or  in  dimensional  form  and  rearranged. 


f(os-V  k)2-cV]  f  v2k2-u2]   *—-   (oj-v  k)2 
I    o         ;km      J    N2      o 


where  the  term  on  the  right  represents  the  coupling  between  the 
two  systems   Now  let  the  perturbation  of  the  elemental  wave 
numbers  be  small: 

k?  =  -~-  +   6,        where  6n  &4 


1   V-c    1         1     V-c 


*      (JO 


k«  -  - —  +  &2  where  62<*  ^\ 

and  keeping  second  order  terms  in  6  we  find: 

3 

62(*2  "  c  V  "  "  ,  .2-7  ,3 

4d  (\/-c) 

If  we  now  let  6*.   — »  0  and  reeall  that  v-c  «  «,  we  see  that: 

1.  The  coupling  is  proportional  to  ratio  of  the  densities 
of  the  two  media  |i,  which  is  a  small  quantity,  and 

2 .  60  is    Hnary  , 

only  if  V  '   c 

tn  the  system  under  consideration 
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Therefore,  c22  -  jk2  and  equations  4.11  reduce  to: 

>x  "  "  jk2a2  *  c23al  (4.12a) 


ki  -*■-    » 


5^-  =  c32a2  ♦  jk^  (4.12b) 


These  equations  can  be  solved  by  assuming  solutions  of  the  form 

X  x 

e    and  setting  the  determinant  of  the  coefficients  equal  to 

zero 


1 
2  2 

*  1,2  =  ±  t(  lc321?  "    22>  }1  S  <  ^  >  (4-13) 


where  k-  -  *f- —  and  k„  =  rr~  ,  the  undisturbed  elemental  wave 
1   v-c       /   v 

numbers . 


Since  they  are  approximately  equal: 

i  4  ■  ^-j^o   where    i  =  (69)  -    — 5 r 

i,Z       Z  Z       4dZ(V-c)Z 

amd  the  coupled  modes  can  be  represented  as  : 
a](x)  -  Axe  lx  +  A2e  2X 

a2(x)    -|^*-     (    (^  +jk1)A1e   lx  +   (  S 2  ♦Jk^e  2X) 
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4.  7 Summary 

The  eigenvalue  problem  has  been  stated,  but  was  not  solved 
because  the  method  of  coupling  of  modes  gave  better  insight  into 
this  problem.   We  applied  coupling  of  modes  to  a  simplified  pro- 
blem and  determined  the  perturbation  of  the  wave  nembers  based 
on  the  effect  of  the  coupling.   In  order  to  properly  apply  the 
boundary  conditions  in  this  case,  the  approach  of  coupling  of 
modes  could  be  followed.   But  it  is  not  clear  that  the  neglected 
waves  can  in  fact,  be  neglected.   Since  at  long  wave  lengths,  the 
boundaries  will  have  a  definite  effect  on  both  modes,  this  must 
be  a ceo. in ted  for  in  the  theory.   The  results  of  Chapter  III  and 
IV  do  provide  a  basis  for  experimentation.   If  we  assume  that 
the  wave  numbers  are  only  slightly  perturbed;  we  can  compute 

from  the  physical  length,  k.   from  equation  4.13  and  cd  and  oi. 

i  r      i 

from  the  dispersion  relationship.  Although  these  values  are  not 
exactly  self-consistent,  they  should  be  approximately  so. 
Experimental  techniques  to  achieve  these  values  are  described  in 
the  next  section. 
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CHAPTER  V 


EXPERIMENTAL  PROCEDURE 


In  order  to  experimentally  test  the  predictions  of  the 
theory  developed,  it  is  necessary  first  to  insure  that  the 
"ingredients"  required  are  available.   The  initial  attempt  at 
assembling  a  device  and  making  test  runs  was  not  successful. 
The  step  by  step  approach  will  be  described  in  detail.   This  con- 
sisted of  several  separate  experiments.   The  first  was  testing 
of  the  plate  structure  to  be  used.   The  second  was  investigation 
of  a  means  of  attaining  supersonic  flow.   Next,  the  two  were  com- 
bined to  test  the  flutter  theory. 

5.1  Testing  of  the  Plate 

There  were  two  experiments  performed  on  the  plate  along.   The 
first  was  an  investigation  of  various  types  of  material  and  vari- 
ous end  conditions   The  test  set  up  is  shown  schematically  in 
Figure  5.1   Due  to  the  magnetic  coupling  between  the  two  cir- 
cuits, the  results  of  this  test  were  qualitative  rather  than 
quantitative.   The  results  are  shown  in  table  5.1.   The  end 
conditions  have  a  noticeable  effect  on  the  response.   The  damp- 
ing or  losses  ffn^condition  but  tightly  clamped  appeared  to  nulify 
any  response 

The  frequency  did  increase  with  tension  which  indicates  that 
the  use  of  membrane  analysis  In  the  coupling  of  modes  section  is 
not  totally  in  error.   The  frequency  of  response  was  reasonably 
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TABLE     5.1 


3  AMPLE 

DIMENSIONS 

EkiO    CONDITIONS 

FReQUeMCY 

RESPONSE 

1     cold   Rolled  iaon 

6*  1 ».030 

CLAPlPEO-CLAMPeO 

SI  (cp«) 
103 

300 

SMALL 
SMALL 
SMALL 

&   HEATTWEATeD  CTKl 

9x|tt.oi? 

CLftMPeO-CLOMP€0 

63 
1  I  Z. 

HeoiOM 

3     SPaiNfii   STGEL 

6«  U-OIS 

CLAMPED-ClAMPEO 

140 

MEDIUM 
SHALL 

SftmCLtt     & 

CLAMPED -PK-i^CO 

no  reivieiofi 

MONE 

S»mLh      2 

CLAMPEO-PIHNED 

7* 
(48 

LM*«E 

ricoiuH 

5$ 


close  to  that  predicted.   (See  example  5.1) 


EXAMPLE  5  1 


where 


2  2  m   ,  n 

m        l       w^ 


P 


D  =  EI,  the  rigidity 

u.  -  mass  per  unit  area 

M  «  longitudinal  harmonic  mode 

length 

a>0  "  650  

2       sec 


*2  *  104  cps 


The  second  experiment  on  the  plate  above  was  to  determine 
if  in  fact  the  second  harmonic  frequency  could  be  excited  when  the 
media  had  the  perdominate  properties  of  a  membrane.   Equipment 
similar  to  that  shown  in  Figure  5.1,  was  used  for  this  part  of 
the  experimentation   To  avoid  the  magnetic  coupling,  one  of  the 
electromagnets  was  replaced  by  a  speaker  for  mechanical  excita- 
tion with  weak  coupling.   The  plate  thickness  was  reduced  to  .001 
inch,  to  more  closely  approximate  a  membrane.   The  equipment  is 
shown  pic tor tally  in  Figure  5. 

In  this  configuration  the  plate  frequencies  were  not  detec- 
table The  membrane  frequencies  gave  excellent  response  and  are 
tabulated  in  Table  5.3   The  variation  between  predicted  and 
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actual  frequencies  can  be  explained.   First,  the  actual  para- 
meters for  the  material  are  difficult  to  determine  exactly. 
Secondly,  at  the  higher  values  of  tension,  there  is  a  possibility 
that  the  grips  or  end  clamps  slipped  slightly. 

The  main  conclusion  to  be  drawn  here  is  that  the  second 
harmonic  can  definitely  be  observed.   The  quality  of  the  response 
decrease  with  increasing  frequency  since  the  higher  modes  are 
more  difficult  to  exert.   The  frequency  expected  when  operating 
the  wind  tunnel  could  be  reasonably  predicted  from  this  type  of 
experiment 

The  responses  at  355  and  720  cps  are  probably  due  to  weak 
coupling  with  the  speaker  or  other  portion  of  the  apparatus,  since 
they  are  independent  of  tension. 

5  2  Experiment  to  Obtain  Supersonic  Flow. 

The  physics  of  supersonic  flow  is  a  very  complicated  sub- 
ject and  is  extensively  covered  in  the  literature,  notably  in 
reference  15.   Only  those  principles  which  are  germane  to  this 
experiment  will  be  touched  on.   In  general,  supersonic  flow  can 
not  be  attained  directly  from  an  air  compressor.   It  is  necessary 
to  modify  the  compressor  output  by  passing  it  through  a  set  of 
nozzles.   For  this  experiment  the  nozzles  were  of  the  convergent- 
divergent  type.   The  design  procedure  is  to  attain  Mach  1  flow  at 
the  throat  of  the  test  section  and  then  expand  to  the  velocity 
desired.   The  calculations  for  this  are  shown  in  example  5.2, 
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based  on  the  set  of  nozzles  shown  in  Figure  5.^ 
EXAMPLE  5.2    (Reference  15,  Table  Bill) 


Inlet  Conditions: 


9     ■-■'   125  r>sia   T  »  70°F 
o  o 


a.   Static  Conditions: 

P 
p   -  ^  -  0.665  -^ 
o        ft 


c  -  49.1  t  -  1158 


o         o        sec. 
Throat  Conditions: 


P*/P     -  0.528 
o 


T*/T     -  0.323 
o 

p*/p     -  0.634 

c*/c     *  0.913 
o 

Test  section  conditions  for  Bach  No.  =2.0 
M*  -  1.633  A/A*  -  1.687         A   t  -  6.0cm2 

P,/P  =0.1275  V  -  1733  ft/sec       .  , 

1   O  **«-!,,-««*.«  1  «si««,* 


P,  =16.1  osia 
Pl/p0  -  0.2300  1 

T./T  »  0  555 
i  o 


d.   Mass  Flow 

A*  4  x60  «  245  cuft/rain. 
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Secondly,  the  flow  in  the  test  section  will  be  greatly 
effected  by  the  build  up  of  the  boundary  layer.   This  has  the 
effect  of  a  convergent  nozzle,  which  tends  to  reduce  to  flow 
velocity  toward  Mach  1.   In  the  case  of  frach  2  flow,  the  velo- 
city will  be  reduced  to  Mach  1  in  a  length  of  31  diameters.   In 
this  experiment,  the  diameter  or  depth  of  the  channel  will  be 
taken  as  one  centimeter  to  allow  sufficient  magnetic  coupling. 
Therefore,  the  length  must  be  less  than  31  centimeters;  an  ini- 
tial length  of  20  centimeters  was  used.   (From  this  length,  the 
thickness  of  the  plate  for  a  desired  natural  frequency  can  be 
computed  as  in  example  5.1).   Based  on  these  calculations,  the 
experiment  shown  schematically  in  Figure  5.3     was  constructed. 

The  pressure  readings  were  taken  with  both  the  plate  in 
position  and  removed.   The  readings  were  taken  with  a  guage  and 
compared  with  a  monometer  for  accuracy.   The  readings  obtained 
are  shown  in  Figure  5.4. 

These  results  are  at  a  variance  with  those  predicted  in 
example  5.2,  but  they  do  indicate  supersonic  flow.   The  variance 
is  due  to  several  factors.   First,  the  ideal  area  ratios  are 
significantly  diminished  by  the  boundary  layer  to  an  "effective" 
area  ratio.   Secondly,  the  capacity  of  the  air  ejectors  was 
insufficient  for  running  the  experiment  as  desired.   Also,  the 
pressure  measurements  were  made  at  the  wall  of  the  channel,  pro- 
bably in  the  boundary  layer,  and  not  in  the  center  of  the  flow. 
The  latter  problem  is  difficult  to  circumvent,  and  in  previous 
experiments  (reference  15,  page  135)  good  agreement  between  actual 
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wall  pressure  and  theoretical  mid-stream  pressures  has  been 
attained.   Using  the  lower  pressures  recorded,  the  speed  of  the 
flow  can  be  approximated.   Based  on  the  readings  at  35  psia  inlet, 
the  flow  is  Mach  1.47. 

5.3  An  Experiment  to  Test  the  Flutter  Theory 

We  have  determined  that  the  ingredients  of  the  theory  are 
attainable.   The  flow  can  be  made  8uPersonici  the  plate  can  be 
excited  in  the  second  mode,  and  the  frequency  of  this  mode  varies 
with  the  tension.  Thus  the  results  of  Chapter  IV  indicating 

that  the  perturbation  on  the  wave  number  k  can  be  estimated  from 

2 

the  coupling  seem  reasonable.   The  rigidity  term  Dk  can  be 

2 
neglected  compared  to  the  membrane  term  M  at  the  frequencies  in 

question.   This  is  especially  true  if  the  very  thin  (.001")  plate 
is  used.   The  first  experiments  with  the  combined  ingredients 
used  a  .015  inch  plate.   The  equipment  consisted  of  the  nozzle, 
a  plate,  a  magnetic  detection  device  and  an  oscilloscope.   The 
presentation  on  the  scope  was  confused,  showing  standing  waves 
and  traveling  waves  and  these  responses  were  at  several  fre- 
quencies.  These  were  probably  excite!  by  the  boundary  layer. 
The  tension  was  varied  and  the  frequencies  of  the  response 
increased  with  increasing  tension.   By  varying  the  pressure  the 
response  would  set  in  at  a  certain  pressure  for  a  given  tension  - 
quantitatively  a  higher  tension  required  a  higher  pressure.  At 
a  given  pressure,  providing  it  was  above  the  threshold,  the  vary- 
ing of  the  tension  was  limited  mechanically.   That  is,  the  tension 


r? 


<3*Ini    i 


-IJ9TO  a;."--;<      flaw 

imtol  a/b  ;  banlaSJa 

1scu  r.c-jq  a  a   ad   rrar   w 


>qxg  nA     f  .g 


©3& 


..aq 


T?J9f    ivarf   sV 
mo  woH   arfT      .aldania 

nooaa  aril   nl   tajloxa 

:en92    s.M'.l    /fj"w 
J -joc,'   a d3   3a/f3 
arrai  vioi39i  -res a  gn  >  arf3 

ia     M  orst:;  ., :  a*'  aquoo  b*»joalaan 


86    «»(.'    flflO    _>J    Bflflflfl 


iaoo  arf:i 

a  bfia  e- 
cwoaTa 
Jt  3  b  asav  at 


.a-xuaaa? 
s;J  ads   ,«1   3 


jvoda   86K 


9   ai   airfT      .oolJaaup 
«K?xa  ariT  auf  ai 

'.  a  baau 
r.-izm  a   .frtaiq  « 
• 

aovaar  gnilavai:!  bna 

•v  a  a 

w  aanoqaai 

sup 

*aa*q  navij  a 

ad*  5o  snJ 


58 


could  be  increased  until  the  plate  failed  in  tension,  and  if  the 
tension  was  decreased  the  plate  would  fail  by  excessive  dis- 
placement.  The  latter  occured  at  low  values  of  tension.   It  was 
felt  that  nothing  of  substance  to  support  the  theory  could  be 
gained  from  further  experimentation  with  this  set  up.   "Flutter" 
did  occur  and  this  did  provide  a  variable  coupling  between  the 
two  magnetic  circuits  and  there  was  an  output  signal. 

However,  this  flutter  could  not  be  attributed  exclusively 
to  the  mechanism  proposed  in  the  theory.   It  was  probably  caused 
by  excitation  from  the  boundary  layer.   There  was  a  component  in 
the  vicinity  of  the  frequency  expected  as  computed  below: 


(1)  length  of  membrane:   k  -  -J  -  .628  in 

(Mfl-l)kvo 

(2)  from  dispersion  equation:   f  -  r =  1080  cps 

a).  ■  495  sec 


3  , 

(3)   from  coupling  of  modes:  t   -  -**« r   =  .25  in 

4d^(v-cr 

In  order  to  prevent  the  end  connection*  from  interferring 
with  the  flow  they  were  removed  from  the  test  section  as  shown 
in  the  figure.   Therefore,  the  flow  is  in  transition  over  a  short 
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section  at  both  ends  of  the  plate.   This  could  also  provide  an 
excitation  which  was  not  predicted  and  would  be  difficult  to 
ana lyze . 

Therefore,  another  experiment  was  designed.   Due  to  a  time 
problem,  it  may  not  be  completed  before  the  due  date  of  this 
thesis.   The  changes  include  using  a  .001  inch  plate  to  more 
closely  approximate  membrane  performance,  refined  end  connectors 
to  insure  a  good  response  and  no  flow  disturbance  by  the  tension 
mechanism.   The  results  of  this  phase  of  the  experimentation  will 
be  reported  in  Appendix  B. 
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CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

6. 1  Conclusions 

The  analysts  of  the  model  does  predict  the  desired  pheno- 
menon.  The  interaction  between  the  flow  and  the  elastic  media 
at  long  wavelengths  will  theoretically  produce  growing  oscilla- 
tions.  The  effects  of  the  longitudinal  boundary  conditions  can 
be  approached  by  the  coupling  of  modes.   In  the  simplified  pro- 
blem considered  in  Chapter  IV,  the  effect  of  the  coupling  on  the 
wave  number  was  shown.   This  approach  could  be  expanded  to  show 
the  influence  of  the  boundaries  on  the  wave  numbers. 

The  experimental  portion  shows  that  the  separate  systems 
can  be  achieved.   The  plate  can  be  excited  in  the  second  mode. 
The  natural  frequencies  of  the  plate  increase  with  tension  and 
the  tension  required  by  the  theory  is  physically  realizable.  A 
supersonic  flow  over  the  plate  can  also  be  attained  in  the  la bora 
tory.   With  the  equipment  used,  the  maximum  velocity  was  1.47 
times  the  speed  of  sound.  The  major  source  of  difficulty 
jicountered  in  combining  the  systems  was  the  excitation  provided 
by  the  boundary  layer.   This  was  not  included  in  the  model  or  in 
the  analysis.   The  excitation  gave  forced  responses,  mainly  in 
the  first  mode,  although  there  was  a  second  mode  in  the  response. 

There  are  several  other  considerations  which  have  not  been 
covered.  With  respect  to  the  elastic  media,  there  will  be  eddy 
currents  induced  when  it  oscillates  in  the  magnetic  field. 
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These  represent  a  loss  In  the  system  which  could  be  significant. 
The  Curie  temperature  will  impose  a  thermal  limitation,  restrict- 
ing efficient  operation  of  the  device  to  relatively  low  tempera- 
tures.  From  the  stand  point  of  strength,  it  is  fortunate  that 
ferrous  materials,  which  are  required  magnetically,  also  show 
good  fatigue  properties. 

The  Tir>st  severe  limitation  on  the  device  efficiency  will  be 
the  losses  incurred  by  the  fluid  in  tr  ansiting  tne  nozzle  and 
interaction  area.   It  was  noted  in  the  experimental  section  that 
the  flow  would  be  supersonic  only  for  a  short  distance.   For 
instance,  the  Mach  1.47  flow  in  the  experiment  would  be  "choked" 
to  Mach  1.0  at  the  end  of  the  test  section  (six  inches),  due  to 
boundary  layer  build  up.  At  the  same  time  the  isentropic  stag- 
nation pressure  would  be  reduced  to  less  than  one  half  its 
initial  values.   This  represents  a  considerable  loss  in  efficiency 
A  more  refined  design  would  improve  on  this  situation,  but  there 
will  be  a  definte  limiting  efficiency  attainable. 

The  subject  of  supersonic  flow  in  short  ducts  is  somewhat 
complex  and  experimental  data  has  not  been  consistent. 
In  general,  the  boundary  layer  thickness  increases  with  lergth. 
But  in  some  cases,  negative  friction  coefficients  have  been 
attained,  indicating  a  decreasing  boundary  layer.   The  effects 
of  the  oscillating  plate  on  the  boundary  layer  have  also  been 
studied  but  there  appears  to  be  no  significant  change  in  boundary 
layer  tl  ickness.  '   ' 
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The  device  itself  may  not  be  the  most  efficient  way  of 
obtaining  energy  from  the  flow.   However,  it  is  possible  that 
to  employ  it  in  tandem     with  a  ^q^   conventional  energy  con- 
verter would  improve  overall  efficiency.  The  important  conclu- 
sion here  is  that  a  mechanism  is  available  for  transfer  of  energy 
between  the  flow  and  the  media.   The  mechanism  itself  is  not  a 
"lossy"  one  -  the  inefficiencies  are  incurred  in  the  scheme 
for  using  the  mechanism  attempted  here.   In  other  words,  this 
paper  has  been  limited  in  scope  to  a  particular  device,  but 
further  research  would  lead  to  other  methods  of  using  this  con- 
cept for  energy  transfer. 

6, 2  Recommendations 

The  study  of  the  longitudinal  boundary  conditions  in  this 
scheme  could  be  expanded  using  the  coupling  of  modes  approach. 
The  coupling  coefficients  considered  here  pertained  only  to  the 
coupling  between  two  waves,  one  in  each  system.   The  coupling 
between  the  other  two  waves,  both  with  the  other  wave  in  its  own 
system  and  the  other  wave  in  the  other  system  were  neglected. 
The  latter  still  appears  reasonable.   But  the  boundary  conditions 
will  provide  coupling  between  the  two  waves  in  one  system.  We 
neglected  this  here  because  at  the  wavelengths  in  question,  the 
rigidity  term  could  be  neglected  and  the  plate  considered  as  a 
membrane.   In  the  general  case  this  would  not  be  so,  and  the  area 
of  longitudinal  boundary  conditions  in  coupled  systems  has  many 
possibilities  for  study. 
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The  combined  experiment  showed  the  ovoaitetanding  effect  of 
the  boundary  layer,  essentially  making  the  desired  response 
if  present.   The  boundary  layer  could  be  reduced  either  by 
thermal  methods  or  by  expanding  the  flow  in  the  depth  dimension. 
If  this  were  doae,  the  correlation  of  experiment  and  theory 
might  be  attained  and  the  mechanism  studied. 

(9  ) 
Ketterer     used  a  fluid  atream  and  a  spring  coupled 

(17) 
electrically  to  study  this  mechanism.   Fraize     studied  a 

device  in  which  a   fin  t   situated  in  a  wind  tunnel  (sub  sonic  ) , 
was  excited  mechanically  to  produce  traveling  waves.   Signifi- 
cant reductions  in  drag  were  recorded.   Replacing  the  mechanical 
excitation  by  magnetic  traction  along  the  length,  using  the  cot^p- 
ling  mechanism  investigated  here,  provides  the  possibility  for 
pumping  of  or  propulsion  through  a  fluid  surrounding  the  fin 

This  would  be  the  inverse  of  the  problem  studied  here,  but 
would  not  have  the  boundary  layer  problem  inherent  in  super 
sonic  flow. 


:  ~.~    a. 


&v»ri 


64 


Appendix  A 
A  COMPARISON  WITH  FLUTTER  THEORY 


Flutter  is  ^defined  as  the  self-excited  oscillations  of  a 

( 18} 
airfoil   in  supersonic  flow.      This  definition  leaves  some 

doubt  as  to  the  physics  involved,  but  since  it  is  restricted  to 

supersonic  flow  it  seems  conceivable  that  the  mechanism  involved 

is  similar  to  the  feedback  situation  described  in  Chapter  III. 

The  problem  of  flutter  is  stated  as  a  partial  differential 
equation,  where  structural,  inertia 1,  and  also  dynamic  operators 
are  combined  to  determine  the  equation  of  motion  for  the  panel. 
Two  approaches  are  followed  from  here  depending  on  the  geometry 
involved.  One  is  to  apply  the  boundary  conditions  to  the  plate 
only.  Since  these  are  two  at  each  end,  a  "fourth  order"  eigen- 
value problem  results. 

The  second  approach  is  to  assume  traveling  wave  solutions , 

eJ     ,  and  assuming  real  values  for  k,  compute  complex  values 

(ya) 
of  a).   Dugyidji     has  used  both  approaches  and  compared  the 

results.    fjgthen  recommends  which  approach  to  use  based  on  the 

geometry  of  the  problem  considered. 

In  terms  of  our  parameters,  the  equation  of  motion  is: 
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D  "t  - p  "t +  2pA    2 +  Gs  k 


2 


o      (M.  wl-2Hg       1_     *i 

ot 


(l.K»)^^X        1-M2      f    *o 
v       a'  a 


where  the  term  on  the  right  is  the  'aerodynamic  pressure'.1  For 

small  M  ,  this  reduces  to 
a 


2    c\?  TSfi 


and  the  equation  is  the  same  as  equation  (2.10)  except  for  the 

structural  damping  term  G  .   In  our  case,  the  aspect  ratio  is 

s 

zero  and  the  first  approach  is  recommended.   That  is:   the  bound- 
ary conditions  on  the  plate  are  applied  and  a  solution  of  the 
form 


t  .  «  ekl*  +  ?2ek2*  +  ?3ek3x  +  *4ek4x  )  e*0* 


is  assumed.   This  forms  the  eigenvalue  problem,  and  the  solutions 
for  cd  ;  co,  are  shown  in  reference  18.  Many  parameters  have  been 
introduced  and  this  leads  to  lengthy  computations  to  compare  the 
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the  results  of  the  two  theories.   For  the  tension,  rigidity,  and 
plate  thickness  employed  in  the  device,  the  solutions  predict  to.  =  100 
and  co  ■  1000.   The  approximation  is  necessary  because  the 
plotted  data  must  be  extropolated  to  include  the  values  of  ten- 
sion specified.   These  results  show  correlation  within  an  order 
of  magnitude  of  those  predicted  in  Chapter  V.   (  to  «=  , 

co,         ) .  Most  of  flutter  investigation  has  been  restricted 
to  determining  the  onset  of  flutter  and  methods  for  avoiding  it»s 
occurence.   Again  using  our  parameters  in  figure  10  of  reference 
18,  flutter  is  predicted  to  occur  when  the  tension  corresponds  to 
making  v  +  ■  v-c.   This  is  surprisingly  consistent  with  the  pre- 
vious results.   The  predicted  flutter  frequency  under  these  condi- 
tions is  on  the  order  of  20to  where  to  equals  the  natural  olate 

o        o  ^ 

frequency  neglecting  tension.   The  previous  results  predict 

©-  ■  10cd  . 
f      o 

There  is  considerable  agreement  between  the  two  results, 

in  spite  of  the  difference  in  boundary  conditions. 
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Appendix  B 


The  previously  described  experiment  was  constructed.   How- 
ever, the  improved  end  conditions  to  lessen  the  flow  disturbance 
were  mechanically  unweildly.   They  did  not  provide  the  firm  con- 
nection desired.   Therefore,  no  further  results  of  significance 
can  be  reported. 
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